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Abstract 

We model the cosmic medium as the mixture of a generahzed Chaplygin gas and a pressureless matter 
component. Within a neo-Newtonian approach (in which, different from standard Newtonian cosmology, 
the pressure enters the homogeneous and isotropic background dynamics) we compute the matter power 
spectrum. The 2dFGRS data are used to discriminate between unified models of the dark sector (a purely 
baryonic matter component of roughly 5 percent of the total energy content and roughly 95 percent gen- 
eralized Chaplygin gas) and different models, for which there is separate dark matter, in addition to that 
accounted for by the generalized Chaplygin gas. Leaving the corresponding density parameters free, we find 
that the unified models are strongly disfavored. On the other hand, using unified model priors, the obser- 
vational data are also well described, in particular for small and large values of the generalized Chaplygin 
gas parameter a. The latter result is in agreement with a recent, more qualitative but fully relativistic, 
perturbation analysis in [l^ . 
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I. INTRODUCTION 



The crossing of different cosmological observations, in particular the anisotropy spectrum of the 
cosmic microwave background radiation (CMBR), the luminosity distance of supernovae of type la, 
gravitational lensing and baryonic acoustic oscillations, indicates that around 95% of the cosmic 
substratum is not directly detectable through electromagnetic emission 

naaa. 

As long as 

one accepts General Relativity (GR) to be valid, the now widely accepted conclusion is that most 
of the substance in the universe must be of non-baryonic origin. This dynamically dominating 
non-baryonic substratum is usually divided into two components: dark matter, a pressureless, 
agglomerating component, being present in local structures like galaxies and clusters of galaxies, 
and smoothly distributed dark energy, an exotic fluid with negative pressure. Dark matter is 
required in order to explain the observed anomalies in the dynamics of galaxies and cluster of 
galaxies, as well as to generate the large scale structures in the universe; dark energy is required in 
order to account for the present stage of accelerated expansion of the universe and for the position 
of the first acoustic peak in the CMBR spectrum. The nature of these dark components remains 
a mystery. For reviews on the subject see 0, 1^. 

Among the host of models that have been proposed for dark matter and dark energy over 
the last years, there are unified models of the dark sector according to which there is just one 
dark component that simultaneously plays the role of dark matter and dark energy. The most 
popular proposal along this line is the Chaplygin gas, an exotic fluid with negative pressure that 
scales as the inverse of the energy density [7|. This ^henomenologically introduced equation of 
state can be given a string theory based motivation [8]. It has also been generalized in different 
phenomenological ways ^ . Another example for a unification scenario for the dark sector is a bulk 
viscous model of the cosmic substratum jlOl]. While the Chaplygin gas model (in its traditional and 



there 



generalized forms) has been very successful in explaining the supernovae type la data 
are claims that it does not pass the tests connected with structure formation because of predicted 
but not observed strong oscillations of the matter power spectrum [2]. It should be mentioned, 
however, that oscillations in the Chaplygin gas component do not necessarily imply corresponding 
oscillations in the observed baryonic power spectrum [l^ . 

The generalized Chaplygin gas is characterized by the equation of state 

P = -— ■ (1) 
For ^ > the pressure p is negative, hence it may induce an accelerated expansion of the universe. 
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The corresponding sound speed is positive as long as a > 0. Recently, a gauge-invariant analysis 
of the baryonic matter power spectrum for generalized Chaplygin gas cosmologies was shown to 
be compatible with the data for parameter values a ~ and a > 3 [l^. This result seems to 
strengthen the role of Chaplygin gas type models as competitive candidates for the dark sector. 
The present work provides a further investigation along these lines. While we shall rediscover 
the mentioned results of [U], albeit in a different framework, we also extend the scope of the 
analysis in the following sense. The authors of [l^ have shown that Chaplygin gas cosmologies are 
consistent with the data from structure formation for certain parameter configurations. Here we 
ask additionally, whether or not the data really favor generalized Chaplygin gases as unified models 
of the dark sector. By leaving the density parameters of the Chaplygin gas and the non-relativistic 
matter component, respectively, free, we allow for a matter fraction that can be different from 
the pure baryonic part. This is equivalent to hypothetically admit the existence of an additional 
dark matter component. In other words, we do not prescribe the unified model from the start. 
Moreover, our study is not restricted to the spatially flat case. The 2dFGRS data are then used 
to test whether or not the unified model, requiring that the matter component describes baryonic 
matter with a density parameter of the order of 5 percent only, is favored. 

Now, a precise estimation of the cosmological parameters using the matter power spectrum is 
very involved, since a detailed discussion of many physical processes like free streaming of neu- 



trinos, electron diffusion, etc. is necessary. We s 
conveniently the BBKS transfer function [17 



1 avoid such more complex analysis by using 
5 Q) 19| to impose the initial conditions. We believe 



that this type of analysis retains the essential features of the process and leads to quantitatively 
relevant results. 

Our study relies on a neo-Newtonian approach which represents a major simplification of the 
problem. In some sense, the neo-Newtonian equations can be seen as the introduction of a first 
order relativistic correction to the usual Newtonian equations 
for cosmology 



2QI, 



21 



22 



211 ]. The neo-Newtonian equations 



23( 1 modify the Newtonian equations in a way that makes the pressure 
dynamically relevant already for the homogeneous and isotropic background. This allows us to 
describe an accelerated expansion of the Universe as the consequence of a sufficiently large effective 
negative pressure in a Newtonian framework. While the neo-Newtonian approach reproduces the 
GR background dynamics exactly, differences occur at the perturbative level. However, the GR 
first-order perturbation dynamics and its neo-Newtonian counterpart coincide exactly in the case 



of a vanishing sound speed 



23l |. One may therefore expect that the neo-Newtonian perturbation 



dynamics reproduces the correct GR results on all perturbation scales at least for small values 
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of the sound speed. For constant equations of state it has been demonstrated that the correct 
large-scale behavior in the synchronous gauge is reproduced 2^. On small scales one expects the 



spatial pressure gradient term to be relevant and the difference to the GR dynamics should be 
of minor importance. Since the observational data correspond to modes that are well inside the 
Hubble radius, the use of a Newtonian type approach seems therefore adequate. 

On this basis our analysis extends previous neo-Newtonian studies to the two-component case. 
One of the components is a generalized Chaplygin gas, the other one represents pressureless matter. 
The advantage of employing a neo-Newtonian approach is a gain in simplicity and transparency. 
While in future work all results will have to be confirmed within GR, we shall ensure already here 
that in the region of overlap between GR and neo-Newtonian dynamics our results coincide with 
the corresponding GR results. Our neo-Newtonian approach reproduces the parameter estimations 
for the unified dark matter/dark energy in [l^ also numerically. Backed up by this success of the 
neo-Newtonian approach we then enlarge the scope of our analysis and test the validity of the 



14l | . Denoting the present value of 



unified model itself by relaxing the unified model priors used in 
the Chaplygin gas density parameter by ^Ico, we admit the total present matter density parameter 
OmO to be the sum of an additional dark matter component with density parameter ^dmo and 
the baryon contribution ri^o, i-e., ^mo = ^dmO + ^feo- Leaving the density parameters free, we 
investigate whether or not the unified model with J7c0 ~ 0.96, J7;,o ~ 0.04 and ^IdmO ~ is favored 
by the large-scale structure data. We mention that a similar investigation using supernova type la 
data reveals that the unification scenario is the most favored one 

Our Chaplygin gas cosmology has four free parameters: the value of a, the present Chaplygin 
gas and dark matter density parameters ilcO a-nd ildmO; respectively, and the present Chaplygin gas 
sound speed Vq. There are two main observational sources concerning the matter power spectrum 



today: the 2dFGRS and the SDSS data sets 



2g|. For reasons to be discussed later, we will 



mainly use the 2dFGRS data. For i^dmo ~ 0, i7{,o ~ 0.04 and OcO ~ 0.96, equivalent to the unified 
model, we obtain a very good fit of the date where very small or very large values of a are preferred. 
This reproduces the GR results of 1^ in a Newtonian context. On the other hand, when OdmOi 
0,co and a are left free, large values of ^IdmO and small values of fico are preferred, thus disfavoring 
the unification model. The same result is obtained when all four parameters (^dm.vi, ^rCi, ct and 
A) are left free. If the curvature is fixed to zero, as indicated by the WMAP results [27], implying 
r^co ~ 1 — ^mo, the predictions do not change substantially and a scenario with almost no dark 
energy is again preferred. In all cases, including those for which the unification scenario is imposed 
from the beginning, the minimum values of the statistical parameter are very similar. This 
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does not seem to allow definite predictions of the model. Any conclusion seems to depend on the 
chosen priors. We compare our results with those obtained from the ACDM model, for which 
the power spectrum test indicates ^dmo ~ 0.25. However, almost no restrictions on the value of 
the cosmological constant are obtained. In fact, the matter power spectrum seems to be a good 
indicator for dark matter but not for dark energy. 

The paper is organized as follows. In section |TI] we recall the generalized Chaplygin gas model 
and the basic equations of standard Newtonian cosmology. In section IIIII we introduce the neo- 
Newtonian framework for the two-component model of a generalized Chaplygin gas and pressureless 
matter and establish the perturbation equations for this system. In section ITVl the power spectrum 
is determined, from which the probability distribution functions for each parameter are obtained. 
Our results are discussed in section [TVl 

II. THE GENERALIZED CHAPLYGIN GAS MODEL 

The generalized Chaplygin gas is characterized by the equation of state ([1]) , implying a negative 
pressure and a positive sound speed as long as j4 > and a > 0. The observational constraints 
from supernova type la data indicate that negative values for a are favored, but the dispersion is 
high enough to allow for a large range of positive values for this parameter [llj]. Negative values for 
a imply an imaginary sound velocity, leading to small scale instabilities at the perturbative level. 
Rigourously, the general situation is more complex: such instabilities for fluids with negative pres- 
sure may disappear if the hydro dynamical approach is replaced by a more fundamental description 
using, e,g., scalar fields. However, this is not true for the Chaplygin gas: even in a fundamental 
approach, using for example the Born-Infeld action, the sound speed may be negative if a < 0. For 
this reason we shall not allow a to be negative. 

The traditional Chaplygin gas model is characterized by a = 1. It is a consequence of the 
Nambu-Goto action parametrized in light-cone coordinates. Through some suitable transforma- 
tions, the light-cone parametrized Nambu-Goto action reduces to the action of a Newtonian fluid 
that obeys the equation of state ([T]) with = 1!^. In this sense, it is somehow natural to con- 
struct a cosmological Chaplygin gas scenario within a Newtonian framework. To be precise, the 
symmetries of the Lagrangian are broken when gravity is included. But this drawback can not 
even be cured by using a relativistic version: in order to preserve the symmetries of the original 
Nambu-Goto action a full string model must be implemented. But the Newtonian approach re- 
mains a reasonable approximation because of the mentioned relation between the Chaplygin gas 
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and the Nambu-Goto action. 

Now, establishing a Newtonian model for a universe in accelerated expansion seems to be 
impossible. In traditional Newtonian cosmology the pressure does not play any role in an isotropic 
and homogeneous universe: the universe evolves always with the scale factor a{t) oc t^/^, implying 
a decelerated expansion. This coincides with the relativistic cosmology for a pressureless fluid. 
The pressure becomes relevant only at perturbative level: there the nature of the fluid is essential 
for the evolution of the density contrast. For the evolution of density perturbations in Newtonian 
cosmology see ref. jisl . The specific application to the Chaplygin gas model has been discussed 
in Q. Let *eteh ..s .ain The New.„„,an eo^olog, is deSned th.o.gh the continuity 

equation, the Euler equation and the Poisson equation [15l |: 

^ + V-{pv)=0, (2) 

V^cj) = AttGp . (4) 

In equations ©-(HI) the pressure appears only in the form of a gradient. Hence, the pressure itself 
does not enter the dynamics of a spatially homogeneous background, i.e., the equations do not 
depend on the nature of the fluid. However, at perturbative level the relevant equation is 

6 + 2-6+l^-A7rGp\6 = , (5) 



a I 



where S = ^ is the density contrast, Sp being a first order fluctuation around the background 
solution, Vg = dp/ dp is the sound velocity and k is the wavenumber of the perturbation. If we 
consider a fluid whose equation of state is given hy p = k/o^, the solution is 



with u = e — ^. 

For the Chaplygin gas model the perturbations initially grow as in the matter dominated uni- 
verse, later they decrease, finally approaching zero, which is the value for the cosmological constant 
model [l^ . 



III. NEO-NEWTONIAN APPROACH 



The drawback of standard Newtonian cosmolo gy, the absence of a pressure term in the back- 
ground dynamics, has been cured in a simple way 2^: in the conservation equation ([2]) one takes 
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into account the work done by the pressure during the expansion of the universe. At the same 
time, the equation for the gravitational potential must be modified in order to render the equations 
compatible. This has been done in references [20|, |21 



22l |. The final equations are 



^ + V-{pv)+pS/-v = , (7) 



dt 



+ = -^-VCP , 8 

ot P + P 



V^(t> = 4ttG{p + 3p) . (9) 

For the case of two non-interacting fluids with energy densities pc and pm and pressures pc and 
Pm = 0, respectively, the equations are: 



^ + V ■{pcVc)+PcV -Vc = , (10) 

^ + ^- .V^- = -^^-V<P, (11) 
dt Pc + Pc 

^ + ^-{p^Um) = 0, (12) 



dt 



+ Vm ■ Vivi = -^4> , (13) 
V^cP = 4ttG{p^ + Pc + 3pc) . (14) 



The subscript m stands for pressureless matter and the subscript c for the (generalized) Chaplygin 
gas component. Considering now an isotropic and homogeneous universe with p = p{t), p = pit) 
and V = ^f, we find 



a \ ^ k SttG 



+72 = -^(Pm + Pc) , (15) 



a J 3 

a inG 



{Pc + Pm + 3pc) • (16) 



a 3 

These equations are identical to the corresponding equations for a homogeneous and isotropic 
universe in GR. In a sense, the neo-Newtonian formulation intends to reproduce the equations of 
GR, but in a Newtonian conceptual framework. 

While there is a complete equivalence between the general relativistic and the neo-Newtonian 
equations in the homogeneous and isotropic background, this is no longer the case at the perturba- 
tive level. As already mentioned, the GR first-order perturbation dynamics and its neo-Newtonian 
counterpart coincide exactly only in the case of a vanishing sound speed [23]- But also for small 
values of the sound speed the neo-Newtonian perturbation dynamics will very likely be a reasonable 
approximation. Since the observational data correspond to modes that are well inside the Hubble 
radius, the use of a Newtonian type approach seems adequate, at least as a first and transparent 
step towards a full relativistic treatment. 
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Defining the fractional density contrasts 



Oc = and dm = — 

Pc Pm 



(17) 



for the Chaplygin gas and matter components, respectively, the first-order perturbation equations 
for the system (fTO|) - p^ are 



6c + {2- 



a l + Uc 



a a 



a a' 



+3« 
a 



+ - 



47rGpc(l + uJc){l + ^vl) \ 6c = 47rGp™(l + uJc)6m (18) 



and 



Sm + 2-6m - 4:-nGpm6m = 4:7TGpm{l + Svl)6c , 

a 



(19) 



where = and = The quantity k'^ denotes the square of the comoving wave vector. 
Dividing the equations (jl8p and (jl9p by Hq and redefining the time as t Hq — > t, these equations 
become dimensionless. In terms of the scale factor a as dynamical variable, the system ()18p -()19p 
takes the form 

Lij^(a) 1 + a 



6'^ + \l+g{a)-- 
[a 1 



+ ujc{a) 



-uJc{a) >6'^ 



9{a) ]_ 
a 



(1 + a)wc(a) + 



1 + a \ , , , auJc{a) k'^ljr 



+l^Ma)[l+u:c{am - 3auJc{a)] }6 
2f[a) 



a \ 1 + uJc{a) 

3 f^mO 



2a3/(a 



a2/(a) 

■[1+Ucia)]6m 



and 



- + g{a) 

a 



6' 



3 n 



mO 



3 



cO 



/i(a)[l — 3aLiJc(a)]5c 



(20) 



(21) 



2a3/(a) 2 /(a) 

where Ih = cHq^ is the present Hubble radius and c is the velocity of light. The prime denotes a 
derivative with respect to a and the definitions 



/(«) 

h{a) 
iVcia) 



o _ ^rno + ^^co[^(a) " 3vl /i^^ja^ 
2a[J7mO + ^coa^h{a) + l^fcoa] 
[1 + (1- A)a-3(^+")]TT^ , 
A 

~/i(a)i+" ' 



with 



aA 



(22) 

(23) 

(24) 
(25) 

(26) 
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have been used. Recall that $7^0 = ^dmo+^fco- For the unified model to be an adequate description 
one expects ^2^0 ~ ^^feo- case the data indicate a substantial fraction of ^dmO: the unified model 
will be disfavored. 



IV. THE POWER SPECTRUM: COMPARING THE THEORY WITH OBSERVATIONS 

The power spectrum is defined by 

V = 5l , (27) 

where Sk is the Fourier transform of dimensionless density contrast 6m- We will constrain the free 
parameters using the quantity 

where V° is the observational value for the power spectrum, is the corresponding theoretical 
result and denotes the error bar. The index i refers to a measurement corresponding to given 
wavenumber. The quantity (j28p qualifies the fitting of the observational data for a given theoretical 
model with specific values of the free parameters. Hence, is a function of the free parameters 
of the model. The probability distribution function is then defined as 

F(x„) = Foe->^'(^")/2 ^ (29) 

where the Xn denote the ensemble of free parameters and Fq is a normalization constant. In order 
to obtain an estimation for a given parameter one has to integrate (marginalize) over all the other 



ones. For a more detailed description of this statistical analysis see reference 
The 2dFGRS 



25[ and the SDSS [26:] are the main surveys to obtain matter power spectrum 
data. The last one covers a larger range of scales but the error bars are more narrow for the former 
one. There are some discussions in the literature concerning the relation between the different data 



28l |. In fact, the use of one or the other or the combination of both may result in different parameter 
estimations. For our model, however, the difference in using one or the other set of data is not 
significant (we have verified this!). Hence, from now on we focus on the 2dFGRS observational 
data for the power spectrum. We use the data that are related with the linear approximation, 
that is, those for which kh^^ < 0.185Mpc^^, where h is defined by Hq = 100 • /ikm/s • Mpc. This 
definition should not be confused with the preceding definition of the function h{a). 

To fix the initial conditions we use the following procedure. The ACDM power spectrum is well 
fitted using the BBKS transfer function [LT^. Then, employing the perturbed equations for the 



ACDM model and integrating back from today to a distant past, say z = 1.000, we fix the shape 
of the transfer function at that moment. The spectrum determined in this way is then used as 
the initial condition for our Chaplygin gas model. This procedure is described in more detail in 
references Q- 

To "gauge" our approach, let us first consider the ACDM model. In the general (non-flat) case 
there are two parameters: f^dmo ^-iid J^ao- In figure [T] we show the two-dimensional probability 
distribution function (PDF) as well the one-dimensional PDFs for the dark matter parameter l^dmo 
and for the cosmological constant parameter r^AO) respectively. From the two dimensional graphic 
it is clear that there is a large degeneracy for the parameter Oaoi while the region of allowable 
values for l^dmo is quite narrow. The degeneracy for the cosmological constant density is less visible 
in the one-dimensional PDF graphic, but it is still considerable. Incidentally, the minimum value 
for the parameter is 0.3822 for O^^o = 0.2387 and Oao = 0.5937, corresponding to an open 
universe. 

The four free parameters to be constrained in our Chaplygin gas model are ildmO) ^cO; ^ and 
a. An analysis with four free parameters is computationally hard. For this reason we shall start 
working with sets of three or two free parameters, fixing the remaining one or two, respectively. 
Only afterwards we consider the most general case in which all parameters are left free. This 
strategy will allow us to check the consistency of the final results. The baryonic component O^q 
is kept fixed in agreement with the nucleosynthesis results. We use the value obtained by the 
recent five years WMAP results, Vly^ = 0.043 (with h = 0.72). We will consider the following 
cases: (i) a spatially flat universe with no separate dark matter component, i.e. ^dmo = 0, a 
baryonic component given by fifoo = 0.043 and a dark sector component OcO = 0.957 - there are 
two free parameters, a and A] (ii) a flat universe with the density parameters free, except for the 
condition ilco = 1 ~ ^dmG ~ with the parameter A fixed {A = 0.15 and A = 0.95); (iii) a flat 
universe with a, A and one density parameter free; (iv) the parameter A fixed (with values 0.15 
and 0.95), while a and the two density parameters are free; (v) all four parameters free. Case (i) is 



the configuration studied in reference [IJ]. Our results for the PDF essentially confirm what was 
obtained in [lj|: the one-dimensional PDF, after marginalizing over A, is higher near a = and 
for Q > 2. For A, the one-dimensional PDF, after marginalizing over a, is initially high, then it 
decreases until A ~ 0.7, subsequently it is increasing again. This behavior is shown in figure [2l 
Considering the two-dimensional distribution, the minimum value for is obtained for a = 3.57 
and ^ = 0, with Xmm ~ 0.378, which is a value slightly better than that for the ACDM model, 
Xmin — 0.382. However, the superluminal sound speed renders this model unphysical. (For a 
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possible modification that could preserve causality, see [ij]). 

In case (ii) we relax the restriction that the pressureless matter component is entirely given 
by baryons. It will turn out that this leads to curious results. For vanishing spatial curvature 
rJcO = 1 ~ ^dmo — ^bo is Valid. As before, ^Ibo is fixed and we fix also A = 0.95. Now, varying 
^dmO) we span a two-dimensional PDF which depends on a and ^^dmo- This two-dimensional PDF 
and the corresponding one-dimensional PDFs for a and JlrfmO) respectively, are shown in figure 
[3l Again, values near a = and for a > 2 are favored. On the other hand, the PDF for il-dmo 
decreases as i^dmo increases. This seems to favor the unification scenario which requires a small 
^dmo- However, if we vary ^IcO instead of f^dmO) we find that the PDF for l^co also decreases as 
flcO increases as is shown in figure HI This seems to lead to the opposite conclusion than in the 
previous case: now the unified scenario which requires a large $7c0 seems to be disfavored. Such 
a contradiction seems to be an artifact of the marginalization procedure as can be seen in the 
corresponding two-dimensional PDFs in figures [3] and U) in the first case the probabilities are high 
near a = and for low values of ^dmo > but at the same time the minimum value for is obtained 
for a = and 0,dmO = 1; on the other hand, in the second case, the probabilities are high near 
a = and QcO = 0, which are also the values for which the minimum of is obtained. We 
conclude that under the given conditions the unification scenario is disfavored. The same results 
are obtained for A = 0.15. In all these cases the minimum value for is essentially the same as 
before. 

To test the previous result, we construct again a three-dimensional parameter space with zero 
spatial curvature but leaving a and A free (case (iii)). Again, the PDF for a is initially decreasing 
but increasing later for a > 2, while the PDF for ^mdo is an increasing function of ^Imdo^ as well 
as the PDF for A increases with A. This result is shown in figure [5l If we now vary fico, its PDF 
is a decreasing function of ^IcO as shown in figure EJ 

The curious fact that the anti-unification scenario seems to be favored is confirmed by another 
three-dimensional parameter space study for which a, ^dmo and 0,cO are left free, while A = 0.95 
(case iv). Notice that the curvature now is arbitrary. As can be seen in figure [71 the PDF for a 
shows the same behavior as in the previous cases, while the PDF for ^IdmO now increases with ^dmo- 
On the other hand, the PDF for f^co decreases with $7c0- The anti-unification scenario is clearly 
favored in this case, and there is no contradiction as for configuration (ii). As in the previous cases, 
the minimum value for is around 0.378. 

Varying all four parameters, all the preceding results are confirmed. The one-dimensional PDFs 
for a. A, i^dmo and OcO are displayed in figure [3 It can be seen that the preferred values are either 
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FIG. 1: The two-dimensional probability distribution function (PDF) for QdmO and Qao (left) and the corresponding 
one- dimensional probability distribution functions for the non-flat ACDM model. In the left panel: the darker the 
color, the smaller the probability. 



PDF 2dFGRS 




FIG. 2: The results for the flat case with flto — 0.043, QdmO ~ and QcO = 0.957, corresponding to the unification 
scenario (case (i)). From left to right: the two-dimensional PDF for a and A (with the same color convention as 
before), the best fitting curve for the power spectrum, the one-dimensional PDFs for a and A. 

a <C 1 or Q > 2, while the probabihty is higher for large values of ildmo and small values of f^co- 

Finally, let us consider the particular case a = 0. In this situation, the neo-Newtonian pertur- 
bation dynamics exactly coincides with that of GR. The ACDM model is recovered for A = 1. The 
results for a = and yl 7^ 1 as well as for a = and A = 1, both for a flat universe, are displayed in 
figure [9l the predictions for ^IdmO a-^s essentially the same as for the ACDM model. For A 1, the 
probability for A grows as A approaches unity. This shows that the method employed is consistent. 
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FIG. 3: The results for case (ii) with flto = 0.043, fldmo = 1 ^ ^co — ^to- From left to right: the two-dimensional 
PDF for a and Q,dmO, the one- dimensional PDFs for a and Q,dmO- 



PDF 




FIG. 4: The results for case (ii) with flto = 0.043, flco = 1 — i^mdo ~ ^to- From left to right: the two-dimensional 
PDF for a and QcO, the one-dimensional PDFs for a and Q^o- 




FIG. 5: The results for case (iii) with ildmo = 1 — ^co — ^^bo- From left to right: the one-dimensional PDFs for a, 
^dmo and A. 
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FIG. 7: The results for case (iv). From left to right: the one-dimonsional PDFs for a, Q,dmO and SlcO- 




FIG. 8: The results for the general case with four free parameters (case (v)). From left to right: the one-dimensional 
PDFs for a, A, Q.cO and f^dmo- 

V. ANALYSIS OF THE RESULTS AND CONCLUSIONS 

In the present work we obtained statistical information about the matter power spectrum by 
comparing the theoretical results for the generalized Chaplygin gas model with the 2dFGRS obser- 
vational data. The free parameters of the model are the equation of states parameters a and A and 
the density parameters ^dmo and f^co- The complete four-dimensional analysis is computationally 
hard but still feasible. We have complemented it by a detailed study of the cases for which only 
two or three parameters are free. This allows us to verify the consistency and the correctness of 
the method employed here. 
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FIG. 9: The results for the flat case with q = 0. From left to right: the one-dimensional PDFs for A, for fldmo with 
A 1 and for QdmO with A = 1. 

If the unification scenario with dark matter and dark energy as a single fluid in a spatially 
flat universe is imposed from the beginning (case (i)), the results of reference 1^ are essentially 
confirmed: there are parameter ranges for which the data are well described by the generalized 
Chaplygin gas model. The probability distribution function for a is high for very small (near zero) 
or very large (greater than 2) values of a. Allowing the parameter A to vary, we find that its 
one-dimensional PDF initially decreases with A, but increases as ^ = 1 is approached. Notice that 
values a > 1 imply a super luminal sound speed and are therefore unphysical (see, however, [l4|). 

A different picture emerges for different priors. Leaving the density parameters for the Chaplygin 
gas and the pressureless matter components free, allows us to test the unified models (pressureless 
matter is entirely baryonic) against models in which there is separate dark matter, not accounted 
for by the Chaplygin gas (cases (ii) and (iii)). We find that the unification scenario is clearly 
disfavored. The PDF is highest in regions with very small values for J7c0 ^-nd large values for ^dmo- 
The behavior for a remains essentially the same as in the previous case. This result is confirmed 
when the condition of a spatially flat universe is relaxed and both density parameters are allowed 
to vary freely (case iv) : the one-dimensional PDFs for a and ilco are decreasing functions of a and 
OcO) respectively, while the PDF for OdmO increases with ^IdmO- Finally, the full four-dimensional 
analysis of the phase space (case(v)) reproduces the results for the lower dimensional cases (ii) - 
(iv). 

What is the origin of these apparently contradictory results? The first aspect to be mentioned 
is that the matter power spectrum data only poorly constrain the dark energy component. Even 
for the ACDM model the matter power spectrum gives information mainly on the dark matter 
component, the dark energy component remaining largely imprecise. It is not by chance that the 
dark energy concept emerged from the supernova data. Our results for the Chaplygin gas model 
show that a large amount of dark matter, different from those described by the Chaplygin gas, is 
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necessary to fit the data. However, the dispersion is quite high. For the flat case with a three- 
dimensional parameter space we find at 2cr, that ^dmO = l^ogi- Another point is the use of the 
neo-Newtonian formalism. However, for small values of the parameter a, the main case of interest 
here, the differences to the full general relativistic treatment are not expected to be substantial. 
Moreover, in the cases of overlap the results of the full theory are reproduced. Finally, possible 
statistical subtleties may influence the outcome of the investigation. But as far as we could test the 
statistical analysis (precision, crossing different information, etc), the results seem to be robust. If 
this is really the case, we must perhaps live with the fact that, while the SNe type la data favor 
a unified model of the dark sector this scenario is disfavored if large scale structure data are 



taken into account, unless specific priors are imposed. 
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